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ABSTRACT 

A generic expectation for gas accreted by high mass haloes is that it is shock heated 
to the virial temperature of the halo. In low mass haloes, or at high redshift, however, 
the gas cooling rate is sufficiently rapid that an accretion shock is unlikely to form. 
Instead, gas can accrete directly into the centre of the halo in a 'cold mode' of ac- 
cretion. Although semi-analytic models have always made a clear distinction between 
hydrostatic and rapid cooling they have not made a distinction between whether or 
not an accretion shock forms. Starting from the well-established Galform code, we 
investigate the effect of explicitly accounting for cold mode accretion using the shock 
stability model of Birnboim & Dckcl. When we modify the code so that there is no 
effective feedback from galaxy formation, wc find that cold mode accretion is the dom- 
inant channel for feeding gas into the galaxies at high redshifts. However, this does 
not translate into a significant difference in the star formation history of the universe 
compared to the previous code. When effective feedback is included in the model, 
we find that the the cold mode is much less apparent because of the presence of gas 
ejected from the galaxy. Thus the inclusion of the additional cold mode physics makes 
little difference to basic results from earlier semi-analytic models which used a simpler 
treatment of gas accretion. For more sophisticated predictions of its consequences, we 
require a better understanding of how the cold mode delivers angular momentum to 
galaxies and how it interacts with outflows. 
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1 INTRODUCTION 

The process of galaxy formation must begin with gas, ini- 
tially distributed rather smoothhr, collapsing to high den- 
sities. Furthermore, to sustain ongoing star formation in 
galaxies requires a continued input of gas over their life- 
times. As such, the question of how galaxies got their gas has 
received a great deal of attention over the history of galaxy 
formation studies. While the initial stages of this collapse 
are purely gravitational (the gas being dragged along by the 
gravitationally dominant dark matter), after halo formation 
hydrodynamic forces come into play and further collapse is 
mitigated by the interplay of gravity, hydrodynamics and 
cooling processes. 

An accretion shock is a generic expectation whenever 
the gas accretes suporsonically as it will do if the halo virial 
temperature exceeds the temperature of the accreting gas 
(Binney 1977). Models of virialization shocks have been pre- 
sented by several authors (Bertschinger 1985; Tozzi & Nor- 
man 2001; Voit ct al. 2003; Book & Benson 2010) with the 
general conclusion that the shock occurs at a radius compa- 
rable to (or perhaps slightly larger than) the virial radius. 
Much debate has occurred over the existence of such shocks 



— their existence has often been assumed in analytic models 
of galcixy formation^ since Rees & Ostrikor (1977). 

Recent work has examined these issues in greater de- 
tail. Motivated by hydrodynamical simulations (Fardal et al. 
2001; see also Keres et al. 2005; Ocvirk et al. 2008; Keres 
et al. 2009), which show that a significant fraction of gas 
in galaxies has never been shock heated, Birnboim & Dekel 
(2003) developed an analytic treatment of virialization shock 
stability. The virial shock relics on the presence of a stable 
atmosphere of post-shock gas to support itself. If cooling 
times in the post-shock gas are sufficiently short, this at- 
mosphere cools and collapses and can no longer support the 
shock. For cosmological halos this implies that shocks can 
only form in halos with mass greater than 10^^ Mq for pri- 
mordial gas (or around IO^^Mq for gas of Solar metallicity). 
These values are found to depend only weakly on redshift 

^ This has always been understood to be an approximation, valid 
only in specific mass regimes: White & Frenk (1991), in discussing 
the fate of the gaseous component of a halo, note that "cooling 
rates may be short enough for the multiphase structure to survive 
the shocks, and it is then unclear how the dynamics of the gas 
component should be modeled." 
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and are in good agreement with the results of hydrodynam- 
ical simulations. 

As a result, in low mass halos gas tends to accrete 
"cold" — never being shock heated to the virial temperature 
and instead raining into the halo as cold clumps along fil- 
aments^. If this gas is to make it into the galaxy it must 
nevertheless lose its energy, either by drag processes in the 
halo or through a shock close to the galaxy which turns its 
kinetic energy into thermal energy which is immediately ra- 
diated away. Halos which do support shocks arc expected 
to contain a quasi-hydrostatic atmosphere of hot gas. The 
structure of this atmosphere is determined by the entropy 
that the gas gains at the accretion shock and that may be 
later modified by radiative cooling (Voit ct al. 2003; Mc- 
Carthy et al. 2007). The transition from cold to hot mode 
accretion is not sharp — halos able to support a shock still 
experience some cold mode accretion (the relative contribu- 
tions of the cold mode decreasing with halo mass). 

The consequences of cold vs. hot accretion for the prop- 
erties of the galaxy forming from such an accretion flow 
have yet to be fully worked out. As Croton et al. (2006) 
have stressed, the absence of an explicit cold mode may 
not be important since the cold gas accretion rate in small 
haloes is limited by the growth of the halo rather than by 
the system's cooling time. In contrast. Brooks et al. (2009) 
demonstrate in hydrodynamical simulations that cold mode 
accretion does allow accreted gas to reach the galaxy more 
rapidly, by virtue of the fact that it does not have to cool 
but instead merely has to free-fall to the centre of the halo 
(starting with a velocity comparable to the virial velocity). 
This results in earlier star formation than if all gas were 
assumed to be initially shock heated to the virial tempera- 
ture of the halo. It is also clear that the situation needs to 
be carefully resissessed in the presence of effective feedback 
schemes that prevent excessive star formation, particularly 
in the high redshift universe. 

In this work we implement a treatment of cold-mode 
accretion into the Galform semi-analytic model of galaxy 
formation, following the methodology of Birnboim & Dekel 
(2003). This will allow us to assess the importance of cold- 
mode accretion for cosmological populations of galaxies 
across a range of redshifts, and to suggest additional studies 
of the cold mode which might improve our understanding of 
its role in the process of galaxy formation. An important as- 
pect is that we are able to explore how feedback and reheat- 
ing of cold gas moderate the cold mode. We note that Cat- 
taneo et al. (2006) have previously explored a simpler imple- 
mentation of cold mode accretion in a semi-analytic model 
of galaxy formation. They parameterized cold mode accre- 
tion by defining a critical mass scale (which had some de- 
pendence on redshift) above which accretion switched from 
cold mode to hot mode. While motivated by the results of 
the Birnboim & Dekel (2003) calculation, this parameteri- 
zation did not capture the full generality of that work as we 
will attempt to do herein. 

The remainder of this paper is arranged as follows. In 
§2 we describe our implementation of cold mode accretion 



^ While this picture seems reasonable on theoretical grounds, it 
as yet has little direct observational support (Steidel et al. 2010). 



in the Galform semi-analytic model. In §3 we present our 
results and, finally, in §4 we discuss their implications. 



2 COLD MODE ACCRETION MODEL 

We implement cold mode accretion in our semi-analytic 
model using the results of Birnboim & Dekel (2003). Birn- 
boim & Dekel (2003) derive a simple criterion for virial shock 
stability which states that 

es < es.crit, where eg = porsA(Ti)/Mo, (1) 

where po is the pre-shock gas density, Ts is the shock radius, 
A(r) is the usual cooling function, wo is the pre-shock in- 
flow velocity of the gas and Ti — (3/16)/imrmj/tto/fe is the 
post-shock gas temperature^. The stability parameter es was 
found to be 0.0126 by Birnboim & Dekel (2003) for a gas 
with adiabatic index 7 = 5/3. At each timestep of our cal- 
culation we compute this criterion for each dark matter halo 
using the prc-virialization model described by Birnboim & 
Dekel (2003) in their §5.1. We assume that a fraction, /hot 
of any gas accreted from the intergalactic medium (IGM) 
during that timestep is added to the usual hot gas reservoir 
of the halo (which is assumed to be heated to the virial tem- 
perature by the virial shock). We adopt a simple model in 
which 

/hot = [1 + cxp ({es.crit - es} /Ae)]-i (2) 

such that /hot — i> 1 for ts <C es.crit and /hot — >■ for es » 
fs.crit with the width of the transition being controlled by the 
parameter Ae. Cooling and infall of this hot mode gas then 
follows the standard treatment for our semi-analytic model 
(see Benson &: Bower 2010). The remaining fraction /coid(= 
1 — /hot) is instead added to a new cold mode reservoir in 
the halo. Gas in this cold mode reservoir is allowed to infall 
onto the central galaxy of the halo at a rate 

-^^infall,cold Tcoldmode A^coldmode/te, (3) 

where Mcoidmode is the current mass of gas in the cold mode 
reservoir, ts is the free-fall time from the halo virial radius 
to the halo centre and Fcoidmodc is a parameter of order unity 
which wc will use to adjust the infall rate. Since cold flow gas 
arrives at the virial radius moving with speed comparable to 
the virial velocity we may reasonably expect Fcoidmode > 1- 
In fact, calculations suggest that cold mode gas reaches the 
center of the halo in around one fifth of the freefall time (Yu- 
val Birnboim, private communication). The metal and an- 
gular momentum content of the cold mode reservoir are also 
tracked and we assume that metals and angular momentum 
are transferred from the cold mode reservoir to the central 
galaxy at rates assuming that the mctallicities and specific 
angular momenta of cold mode reservoir gcis elements are all 
equal. 

If a halo merges with a larger halo, its hot and cold 
mode reservoirs are assigned to the hot and cold reservoirs 

^ The left hand side of this expression is equivalent, to order 
of magnitude, to the ratio of sound-crossing and cooling times 
in the post-shock gas. This provides some physical insight into 
this condition: if the post-shock gas can cool too quickly sound 
waves cannot communicate across the halo and thereby form a 
hydrostatic atmosphere which can support a shock front). 
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of the larger halo in just the same way, i.e. the gas from 
the infalling halo is treated just as if it were accreted from 
the IGM. (In practice this process happens gradually as ram 
pressure forces remove gas from the satellite halo. However, 
in most cases the transfer of gas from satellite to host hap- 
pens rapidly after accretion of the satellite.) Our model in- 
cludes outflows of gas from galaxies, driven by supernovae 
(SNe) explosions. Any such reheated gas is always added to 
the hot mode reservoir (after a delay of order a halo dynam- 
ical time sec Benson & Bower 2010). The rest of our semi- 
analytic model remains unchanged and follows the methods 
described by Benson & Bower (2010). 



3 RESULTS 

We find that values of ts.crit ~ 0.0126 (the value found by 
Birnboim & Dokel 2003), Ae = 0.01 and rcoidmodo ~ 2.5 give 
results which are in rcEisonable agreement with those found 
in numerical simulations. It should be noted that these corre- 
spond to a very simple model with little fine-tuning — es,crit is 
as predicted by a simple analytic argument, Fcoidmodc = 2.5 
implies that the cold reservoir is depleted on a fraction of 
the freefall timescale as expected and Ae = 0.01 implies a 
rapid transition from hot to cold mode. 

In particular, Fig. 1 shows a comparison with the results 
from hydrodynamical simulations as reported by Kcros ct al. 
(2009). The fraction of halo gas currently in the cold mode 
reservoir is indicated by the black points as a function of halo 
meiss and for four different redshifts. The blue lines indicate 
the median cold mode fraction at each halo mass, while the 
red lines indicate the corresponding hot mode fraction. For 
comparison, the magenta lines show the cold mode fraction 
as reported by Keres et al. (2009). For reference, the green 
lines indicate the mass at which halos transition from rapid 
to slow cooling (a distinction which has been present in semi- 
analytic models since their inception). Clearly this occurs at 
a lower mass scale than the critical mass for shock formation. 

In the left-hand panel we show results from our model 
with SNe-driven outflows switched off. The simulations of 
Keres ct al. (2009) do not create strong outflows which may 
be expected to significantly alter the nature of cold vs. hot 
mode accretion. Therefore, the fairest comparison to their 
results is made with outflows removed from our model. Our 
model produces results in reasonable agreement with those 
of Keres et al. (2009), correctly finding the characteristic 
transition mass as a function of redshift and producing ap- 
proximately the correct width of the transition at « = 0. 
Our model transitions somewhat too quickly at higher red- 
shifts, which may be indicative of a break-down in the as- 
sumption of spherical symmetry in the model of Birnboim 
& Dekel (2003), an assumption which should be more valid 
at low redshifts. In the right-hand panel of Fig. 1 we switch 
SNc-drivcn outfiows back on in our model, under the as- 
sumption that gas driven out of galaxies will populate the 
hot component of the halo. Wc find that the cold mode frac- 
tion is greatly suppressed in this case, since the halos now 
contain significant quantities of reheated gas. This assump- 
tion could be incorrect, with outfiowing gas condensing into 
small, dense clouds; our goal here is merely to indicate how 
important outflows could be in establishing the cold/hot ac- 
cretion balance in halos. 



In Fig. 2 we explore the consequences of cold mode ac- 
cretion for the build up of the galaxy population. The left- 
hand panel shows the volume averaged star formation rate 
in models with and without cold mode accretion (blue and 
red lines respectively) in cases with and without SNe-driven 
outflows (thick and thin lines respectively). Of course, the 
models lacking SNe-driven outflows produce far too high star 
formation rates as expected, but are shown for reference. 
Comparing the blue and red lines, the models have almost 
identical star formation rates at high redshifts irrespective 
of the inclusion of cold flows both when SNe-driven outflows 
are included and when they are not. At low redshifts, the 
models including cold mode accretion have higher star for- 
mation rates. Wc note that the model without cold mode ac- 
cretion and with SNe-driven outflows has been constrained 
to flt a variety of observational data, including those shown 
in this Figure. The addition of cold mode accretion worsens 
the agreement with the data, but this should not be taken 
as evidence against cold mode accretion — it is highly proba- 
ble that adjustments in other parameters could restore good 
agreement in models with cold mode accretion. Comparing 
thick and thin lines, the flgure emphasizes the importance 
of using the correct galaxy formation physics to assess the 
importance of cold flows. The effects of stellar winds and 
outflows are far more important (but less certain) than the 
distinction between cold flows and rapid cooling. 

The right-hand panel of Fig. 2 shows the fraction of 
baryons (averaged over the entire universe) in the form of 
stars (dashed lines) or stars and interstellar medium (ISM) 
gas (solid lines) as a function of redshift for the same set 
of four models. Considering first models with no SNe-driven 
outflows, the fraction of baryons in stars is identical in mod- 
els with and without cold mode accretion at high redshifts 
as expected from the star formation rate. However, the frac- 
tion of baryons in stars plus ISM gas is actually higher in 
the model including cold mode accretion. This occurs be- 
cause cold mode accretion is more efficient at getting gas 
into the galaxy phase but, in our model, results in galax- 
ies with longer star formation timescales due to a higher 
angular momentum content (and consequently longer disk 
dynamical times). This is illustrated in Fig. 3, where we 
show the dynamical times of galactic disks as a function of 
halo mass at z = 6. Dashed lines compare the effect of the 
cold mode on dynamical time in the absence of SNe-driven 
outflows, while solid lines are for models that include SNe- 
driven outflows. In high mass haloes, the cold mode increases 
the disk dynamical time by a factor of three or more. In 
models that include SNe-driven outflows the effect is weaker 
and the difference in total gas content is only signiflcant at 
higher redshifts. By 2 = 0, in models with SNe-driven out- 
flows, the total mass of baryons in stars and ISM gas differs 
only slightly between models with and without cold mode 
accretion. 



4 DISCUSSION AND CONCLUSIONS 

We have described a simple implementation of cold-mode 
accretion in a semi-analytic model of galaxy formation, us- 
ing a previously proposed analytical model with no modi- 
flcations. Without any adjustment this model provides an 
excellent match to results from numerical simulations — this 
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Figure 1. Tlie fraction of diffuse (i.e. non-galactic) gas in halos in the cold mode is shown by the black points as a function of halo mass 
and for four different redshifts as indicated in each panel. The solid blue line shows the median cold mode fraction at each halo mass. 
The red line shows the corresponding median hot mode fraction. The magenta lines show the equivalent cold mode fraction from the 
simulations of Keres et al. (2009). Green lines indicate the transition from rapid to slow cooling in models without cold-mode accretion 
(halos in the rapid cooling regime are assigned a cold fraction of 1, those in the slow cooling regime are assigned a value of 0). Left-hand 
panel: Results when SNe-driven outflows are not included. Right-hand panel: Results when SNe-driven outflows are included. 




Figure 2. Left-hand panel: The volume averaged star formation history as a function of redshift. Green points show observational 
estimates from a variety of sources as compiled by Hopkins (2004) while magenta points show the star formation rate inferred from 
gamma ray bursts by Kistler et al. (2009). Red lines show results for models with no cold- mode accretion, while blue lines are for models 
including cold-mode accretion. Thick lines indicate models including SNe-driven outflows, while thin lines are for models without such 
outflows. Right-hand panel: The fraction of baryons locked up into stars (dashed lines) and stars plus galactic gas (solid lines) as a 
function of redshift. Thick and thin lines correspond to models with and without SNe-driven outflows. Red lines show models without 
cold-mode accretion, while blue lines correspond to models that include cold-mode accretion. 



Cold Mode Accretion in Galaxy Formation 5 




10'" 10" lO's 



Figure 3. The median dynamical time of galactic disks as a func- 
tion of halo virial mass at z = 6. Dashed lines show models with 
no SNe-driven outflows while solid lines include SNe-driven out- 
flows. Red lines show results for models with no cold-mode accre- 
tion, while blue lines are for models including cold-mode accre- 
tion. 



could no doubt be improved with some fine-tuning of the 
model, but our aim here was to demonstrate that the simu- 
lation results can be encapsulated by an easy to implement 
model. 

We find that the inclusion of SNe-driven outfiows has a 
dramatic efi'ect on the fraction of cold mode gas present in 
dark matter halos — much more so than switching the cold 
mode accretion on or off — suggesting that simulations must 
account for such outflows before they can make robust pre- 
dictions for the properties of cold mode gas. The inclusion 
of cold mode accretion makes little difference to luminosity 
functions and galaxy sizes at 2; = 0, implying that results 
from earlier semi-analytic models which did not treat cold 
mode accretion are still valid. Without SNe-driven outflows 
the cold mode results in a significant increase in the mass 
and luminosity of brighter galaxies which are able to gain 
some mass through the cold mode even when their hot mode 
has been effectively shut down by feedback from active galac- 
tic nuclei (AGN) . 

Cold mode accretion is found to be effective at getting 
gas into the galaxy phase at high redshifts, potentially al- 
lowing for high rates of star formation during these epochs. 
However, we find that the star formation rate at z ^3 is 
mostly unaffected when we include cold mode accretion due 
to the fact that the galaxies that form are larger and lower 
density than they would be if cold mode accretion were ne- 
glected. We caution that the treatment of angular momen- 
tum delivery to galaxies via the cold mode is therefore of 
crucial importance to assessing its impact on star forma- 
tion rates. To date, this has not been studied in detail in 
hydrodynamical simulations but should be in order to re- 
flne our understanding of how the cold mode affects galaxy 
formation. 

In summary, cold mode accretion should be explicitly 



accounted for in semi-analytic models — and will require 
some retuning of parameters to restore good flts to observa- 
tional data — but does not seem to qualitatively change our 
picture of galaxy formation, at least at the coarse-grained 
level studied here. More detailed results (luminosity func- 
tion shapes, distribution of sizes and formation times) will 
require a signiflcantly more detailed implementation of the 
cold mode, including how it delivers angular momentum to 
galaxies (e.g. Navarro et al. 2009). Calibration from numeri- 
cal simulations should improve the accuracy of cold mode re- 
sults. However, we have seen that the inclusion of cold-mode 
accretion in to our semi-analytic model does not produce any 
surprises. The changes in star formation rates and stellar 
mass fractions, particularly in the low redshift universe are 
reassuringly small, and much smaller than the differences in 
galaxy properties created by changes to the feedback effects 
of galactic winds and AGN. 
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